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To assure a real penetration of the enormous intermittent and dispersed naturally renewable resources 
especially solar and wind, optimal sizing of hybrid renewable power generation systems provide evidence 
to be indispensable. A determinist technique has been developed to optimize the annual capital cost and 
the levelized cost of energy of a stand-alone hybrid generation system. Generation and storage units for 
three systems (wind, photovoltaic, and hybrid photovoltaic-wind turbine) are optimally sized in order to 
meet the annual load and minimize the capital annualized cost. A case study is conducted to analyze one 
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1. Introduction 


One of the major challenges confronting users and designers of 
wind and solar energy systems is the casual, unpredictable nature 
of the energy sources |1]. In one hand, the photovoltaic (PV) stand- 
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alone system is an expensive option and it depends on the 
variation of sunshine periods. In the other hand, the stand-alone 
wind energy conversion system can provide an utilizable energy in 
a portion of the time during the year due to high cut-in wind 
speeds which range from 3.5 to 4 m/s [2] and due to the significant 
fluctuations in the magnitude of wind speeds from 1 h to another 
throughout the year. Hence, a stand-alone solar energy or wind 
energy system properly cannot match the electricity demand of 
consumers due to the high monthly and seasonal variations and 
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the unattended operation mode for extended periods of time. For 
example, in the winter seasons, the solar radiation is at its lowest 
level while the wind speed is usually at its highest level. Moreover, 
during the nights, solar energy cannot be utilized while the wind 
energy may be exploited efficiently. Hence, simultaneous utiliza- 
tion of multiple energy resources by combining solar and wind 
energy conversion systems enhances greatly the certainty of 
meeting load demands. The combined systems, called also hybrid 
of energy systems (HES), are becoming more and more attractive 
and being extensively used for substitution of fossil-produced 
energy, and eventually to decrease air pollution [2]. 

In reality, the electrical demand can be random that’s why 
stand-alone HES may suffer from intermittent breakdown, which 
impacts its supply reliability. This may make the HES unreliable 
relatively compared to the traditional supplies of electrical energy. 
Therefore, the use of some energy storage devices seems so 
necessary to provide a high reliability and avoid gross overdesign 
of the solar and wind system [3,4]. Another advantage of the HES 
is the possibility to reduce the required energy storage capacity 
when one of the optimum combinations of photovoltaic and wind 
energy is used [5,6]. These devices store the excess of electricity 
demand and subsequently meet the load demand in shortage of 
HES. The conventional lead-acid battery is the most general energy 
storage device at the actual time [7]. Generally storage systems are 
very expensive and need to be reduced to their minimum possible 
storage size for a renewable energy system in order to be cost 
effective [8]. However, a hybrid solar-wind power system 
improves the overall energy output and reduces energy storage 
requirements [8]. 

The output generation of the PV and wind turbines depends on 
the wind speed and solar radiation variations, that’s why the 
batteries are used to store the surplus energy produced and 
discharge it when the renewable sources are not available [9]. 
The converter operates at maximum power point tracking mode to 
extract the maximum possible energy, regardless of the variation 
in weather conditions. Various operating strategies can be used, 
depending on the wind and solar availability and the controller 
topology. For instance, when solar power is not available at night, 
wind and battery backup will be used to supply the load demand. 
Moreover, when the output power generated by the two sources 
and the stored power in the battery bank are insufficient to meet 
the load power consumption, DC backup generators can be used to 
increase the reliability of the system. 

For any given load condition, an accurate sizing of PV, wind 
and battery lead to a successful operation of a hybrid energy 
system [10]. 

Hybrid power systems are mainly used in remote and rural 
locations where the cost of supplying power is high. Besides, the 
daily load profiles are an important parameter in the hybrid system 
design, as it characterizes the power demand that must be met by 
the available renewable sources and backups. As the renewable 
sources vary throughout the day, the excess power is stored and the 
batteries will act as auxiliary supply to meet the load demand. 

Nevertheless, the hybrid energy system presents some pro- 
blems due to its increased complexity in comparison with single 
energy systems. This complexity makes the evaluation of the 
hybrid energy systems more difficult. The solar radiation and 
wind speed being highly location dependent, the sizing of such 
hybrid systems requires comprehensive analysis of these variables 
for a given site in relation to the system cost for different 
combinations. Hence, it is very important to determine the levels 
of the energy resources at which the response achieves its 
optimum. The optimum design parameters depend on the objec- 
tive function which could be either at its minimum or a maximum 
values. In this study, we aim to develop an interesting optimiza- 
tion method for a PV/wind hybrid energy system. The objective 


function has to depend on PV and wind turbine power output and 
battery capacity. 


2. Inspiration and background 


In the past, a several surveys [11,12] were based upon a 
particular scenario with a definite set of design values yielding 
the optimum solution. Unfortunately, a lot of last approaches 
provided the optimum solutions without the aptitude to supply 
a general appreciation about how the total system cost changes as 
a function of the parameters design size. Recently, numerous 
research groups have carried out the optimization of autonomous 
hybrid energy systems. Borowy and Salameh [13] proposed an 
algorithm to optimize a photovoltaic-array with a battery bank for 
a standalone hybrid PV/wind system basing on a long-term hourly 
solar irradiance and peak load demand data. Later, Borowy and 
Salameh [14] optimized a PV/wind system combined with a 
battery bank taking into account the cost of the PV modules and 
battery units. Markvart [15] constructed a graphic technique to 
optimize the size of the PV/wind energy system basing on the 
monthly average solar and wind energy values of the south of 
England. Bagul et al. [16] used the long term data of wind speed, 
irradiance and ambient temperature measured every hour during 
30 years and the load specifications for a typical New England 
house. They proposed a technique to determine the optimum 
relationship between the number of PV panels and the number of 
required storage batteries for the stand-alone hybrid wind-photo- 
voltaic system, to match a certain loss of power probability. 
Morgan et al. [17] enhanced the sizing and the optimization of 
the autonomous renewable energy systems by predicting the 
battery state of voltage (SOV) rather than its state of charge 
(SOC). In fact, the simulated algorithm permitted them to forecast 
the hybrid energy system performance as a function of battery 
temperature. Celik [18] used the long measured hourly weather 
data of eight years in order to make a techno-economic investiga- 
tion and optimization of a PV/wind hybrid energy system basing 
on 1996 weather data from TyB site of Cardiff, UK. The perfor- 
mance of the HES was evaluated in terms of solar and wind energy 
fractions, battery storage capacity and system cost. Yang et al. [19] 
utilized an optimization approach basing on the loss of power 
supply probability (LPSP) model for a PV/wind hybrid system in 
order to assess its reliability. Mitchell et al.[1] used a heuristic 
model, based on seasonal averages of wind, solar and load 
consumption in order to optimize HES for both grid-connected 
and stand-alone applications. Their model was able to minimize 
capacity of batteries and stand-by energy import. Tina and al. 
developed probabilistic methods which incorporated the irregular 
nature of the resources and the electrical demand in order to 
eliminate the requirement for time-series data and evaluate the 
long-term performance of HES [20]. Ashok [21] suggested a model 
to minimize the life cycle cost of a hybrid energy system and 
proposed a general model to find an optimal combination of 
energy components for a usual rural society. Ekren [22] introduced 
the response surface methodology as a numerical tool to the size 
optimization of PV/wind hybrid energy system for a given time- 
varying hourly load demand that is assumed to be the electricity 
consumption of the global mobile communications station at Izmir 
Institute of Technology Campus Area, Urla, Turkey. The achieved 
simulation was based on the hourly mean solar radiation and wind 
speed data for the period 2001-2003 registered at the meteor- 
ological station of Izmir institute. The Response Surface Methodol- 
ogy, used by Ekren, led to a better understanding of the accurate 
relationship between input variables such as PV size, wind turbine 
rotor swept area and battery capacity, and output variables such as 
the hybrid energy system cost. 
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Fig. 1. Typical design of hybrid energy system. 


3. Hybrid energy system 


The hybrid system relies on solar and wind energies as the primary 
power resources, and it is backed up by the batteries (see Fig. 1). 

Here, batteries are used due to the stochastic characteristics of the 
system inputs, solar radiation, wind speed and the electricity con- 
sumption of the load. The potential energy resources such as solar and 
wind are not controllable variables and the behavior of these variables 
are non-deterministic. Therefore, probability distributions are specified 
in order to carry out a random-input simulation. Hourly solar radia- 
tion’s and wind speed’s theoretical distributions are fitted to a suitable 
used distribution. In the simulation model, hourly data are used. One 
of the hybrid system’s simulation model assumptions is that the input 
variables do not change throughout an hour. The length of each 
simulation run is considered as 20 years of the economic life which 
consists of 365 days/year, 24 h/day, in total 175,200 h. For each run, 
five independent replications are completed. In the simulation model, 
the annualized capital cost C; is used to compare two or more 
alternative configurations since it is a popular and useful economical 
technique [23]. 


3.1. Solar radiation 


The climatic conditions of the Monastir city are favorable for 
operation in various renewable applications with an annual sunshine 
hours and an average daily solar radiation received on horizontal plane 
equal to 2775 h and 4767 kW h/m?, respectively. We note that the 
sunshine duration accounts the number of hours for which the 
measured solar irradiance goes over 120 W/m2. The sunshine duration 
can be expressed other wisely in term of sun cover rate. 


Taher et al. [24] recommended that the optimal yearly tilt angle 
of a fixed solar collector, faced due to the south, in the Monastir 
city is close to 0.9 times its latitude (33°). Basing on this result, in 
this study we will use the average hourly solar radiation received 
on a tilted plane with an inclination of 33° (Fig. 2). 


3.2. Wind speed 


Monastir city is known by a daily average wind speeds between 
4.7 and 7.2 m/s and both stable and high temperatures during the 
hot and the sunny days. The average hourly wind speed variation 
gives information about the availability of appropriate winds during 
the entire 24 h of the day. The hourly wind speed profile which was 
recorded during the period 1983-2006 is shown in Fig. 3. 


3.3. Hybrid energy system modeling 


3.3.1. PV model 

Knowledge of global solar radiation, using either measured 
values or estimated ones basing on reliable models, is necessary to 
PV performance modeling. 


3.3.1.1. Global solar radiation 

3.3.1.1.1. Beam radiation. The beam radiation on tilted surface, Gp, 
can be estimated by multiplying its value received on a horizontal 
surface, Gp, by a geometric factor R, which depends on the zenith and 
incidence angles. R,is given by the following the expression [25-27]: 
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Fig. 2. Average hourly solar radiation of Monastir city received on a tilted solar collector. 
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Fig. 3. Long term average wind speed of Monastir city (1983-2006). 


Thus, the beam radiation on tilted surface Gp is expressed 


as: 


Gre = Gp x Rp = Gp x 


cos 0 


cos 0, 


(2) 


3.3.1.1.2. Diffuse radiation. The diffuse radiation component is 
the only difference among the models appears in the assessment 
of the sky-diffuse component [28-33]. 

It has been pointed out that the sky diffuse component is 
considered to be the largest potential source of computational 
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error [34,35]. In fact, the most complex problem of estimating the 
diffuse irradiance on a tilted surface, G,,, is related to the fact that 
it does not have an isotropic distribution in the sky which is not 
uniform over the time. Many authors investigated the diffuse 
radiation component basing on various measurements from sev- 
eral stations. In order to evaluate this component, they used 
various models classified basically on three types: isotropic, 
circum-solar and anisotropic. 

3.3.1.1.3. Isotropic diffuse radiation. The isotropic models pre- 
sume that the intensity of diffuse radiation is uniform over the 
entire hemisphere sky dome which could be true in the case of 
overcast skies. Hence, the diffuse radiation incident on a tilted 
surface depends on the fraction of the sky dome seen by it [36]. 

Dorta [37] revealed that many models calculating solar radia- 
tion on an tilted surface with any particular orientation have been 
performed and the isotropic diffuse sky model, Hottel-Woertz- 
Liu-Jordan [38,39], is the most used. According to this model, the 
diffuse component can be given as follows: 


Gat = Ga x (SE = A 


(3) 
where £ is the slope of the considered surface. 

However, despite the popularity of the isotropic models, 
theoretical as well as experimental results have shown that the 
simplified calculations of isotropic model do not give a real picture 
of the insolation conditions [37]. In fact, they showed the poorest 
performance and should not be used for estimating the diffuse 
radiation on a tilted surface [40]. 

3.3.1.1.4. Circum-solar diffuse radiation. The circum-solar model 
presumes that all the radiation comes from the direction of the 
sun and its surroundings, and, therefore, the diffuse component is 
calculated analogously to the direct component [25]: 


Gat = Ga x Rp (4) 


Unfortunately, this model could only be applied in the case of 
completely clear skies, generally overestimates the diffuse compo- 
nent of the radiation [26]. 

3.3.1.1.5. Anisotropic diffuse radiation. The anisotropic model 
assume the anisotropy of the diffuse sky radiation in the circumsolar 
region (sky near the solar disc) and the isotropy diffuse component 
distribution in the rest of the sky dome. This model is widely accepted, 
it is used both for clear and cloudy or partly cloudy days [26]. Based on 
the above, an important number of models considering anisotropic 
distribution of the diffuse irradiance in the sky are proposed [41-45]. 
In this context, one of the most confirmed anisotropic sky model is the 
HDKR, Hay-Davies-Klucher-Reindl, model [37]. This model assumes 
that there are three components of the diffuse solar radiation: an 
isotropic component which comes from all parts of the sky equally, a 
circumsolar diffuse component, which is concentrated in the sky near 
the sun, and a horizon brightening component which emanates from 
the sky near the horizon (Fig. 4). 

According to The HDKR model, the diffuse component radiation 
incident on a tilted surface is as follows: 


Gar = Ga(RpAi) +1 Ay a A) h sin? (5)] = 


where f is the factor which accounts the fact that more diffuse 
radiation comes from the horizon than from the rest of the sky, 
this modulating correction factor of diffuse radiation includes 
influence of cloudiness and it is expressed by: 


Gp 
=\ = 6 
f 7 (6) 


Ai is the anisotropy index used to estimate the amount of 
circumsolar diffuse radiation, also called forward scattered 
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Fig. 4. Schematic view of the solar radiation distribution on a tilted. surface 
according to anisotropic diffuse sky model. 


radiation. This factor is expressed as a ratio of the beam radiation 
on a horizontal ground surface to the extraterrestrial radiation. Its 
high value enhances the contribution of circumsolar diffuse 
radiation. The anisotropy index is given by: 


_G, 


Aj 
Go 


(7) 


3.3.1.1.6. Ground reflected radiation. All models assume that the 
ground reflected component is isotropic and it is expressed as 
follows [46,47]: 


Ga one =A 


(8) 
where p is the albedo coefficient of the ground. 

3.3.1.1.7. Global radiation. The global radiation on a tilted surface, 
Gr, is the sum of the diffuse radiation, Ga , the beam radiation G+, and 
the ground reflected radiation, G,. Therefore, the incident global radi- 
ation on tilted surface is given by the following expression [48-53]: 


Gt = Gbr + Gat + Gr (9) 


In this study, the HDKR model will be applied assuming that the 
global radiation on a tilted surface is given by the following formula: 


Ge = (Gp + GgAi)Rp +Ga(1 -ARa i +f sin? (5)] +Gp,R, (10) 


3.3.1.2. PV generation. The PV generation is mainly dependent on 
the global incident radiation (which was overviewed in the 
previous section) and the PV cell temperature. 

3.3.1.2.1. PV cell temperature. It is well known that the PV cell 
temperature has an important effect on the PV output power. This 
temperature can be the same as the ambient temperature during the 
night but it can go over the ambient temperature by 30° or more in 
full sun. Hence, it is necessary to calculate the PV cell temperature in 
order to account for this meaning effect. For this reason, we have to 
begin by establishing the energy balance between, on one hand, 
the solar energy absorbed by the PV array, and on the other hand, the 
electrical output plus the heat transfer to the surroundings. 
The balance on unit area of the module, which is cooled by losses to 
the surroundings, can be written as [25]: 


TAGT =Nc Gr +U,(Tc— Ta) (11) 
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where 7, @, Gr, Nc, Ur and Ta are respectively the solar transmittance 
of the PV array, the solar absorptance of the PV array, the global solar 
radiation striking the PV module, the electrical conversion efficiency of 
the PV array, the coefficient of heat transfer to the surroundings and 
the ambient temperature. 

According to The above equation, the PV cell temperature can 
be expressed as follows: 


ee ee or(F) (1 tc) (12) 


Ta 


To estimate the value of (ta/U,), we report the nominal operating 
cell temperature (NOCT), which is defined as the cell temperature that 
results at an incident radiation of 0.8 kW/m/?, an ambient tempera- 
ture of 20 °C, no load operation (meaning 7-=0) and an average wind 
speed equal to 1 m/s [54,55]. We can substitute these values into the 
above equation and solve it for (tæ /U;): 


ta _ Tcnocr— Tanocr 
: : 13 
UL Gr.ocr a3) 


where Tcnocr, Tanocr and Gryocr are respectively the nominal 
operating cell temperature, the ambient temperature and the solar 
radiation at which the NOCT is defined. 

If we presume that (ta/U;) is constant, we can substitute this 
equation into the cell temperature equation as follows: 


(Tcnocr — Tanocr) Nc 
fe= Tik Gr( oe ) (1 a) (14) 


We suppose that the PV array always operates at its maximum 
power point, which means that the cell efficiency is always equal 
to the maximum power point efficiency mp: 


Nc = mp a5) 


In this order, we can substitute nc by 77, in the relation (14) 
which can be written therefore as follows [56]: 


G 
Tc = Ta+(Tcnocr— Tanoct) (z z) (1 = Tee) (16) 


However, mp depends on the cell temperature Tc, we can 
assume that the efficiency varies linearly with temperature 
according to the following equation: 


"mp =Nmpstcl ap(Tc Tcstc)] (17) 


where mp secs Ap and Teste are respectively the maximum power 
point efficiency under standard test conditions (a radiation of 
1 kW/m, a cell temperature of 25 °C, and no wind), the tempera- 
ture coefficient of power and the cell temperature under standard 
test conditions. We note here that the temperature coefficient of 
power (ap) is negative which means that the efficiency of the PV 
array decreases with increasing cell temperature. 

Finally, by substituting this efficiency Eq. (17) into the above 
cell temperature Eq. (12), the cell temperature becomes: 


Ta+(Tcnoct — Tanocr)(Gs/Gs,nocr)[1 — (mp ste/TOU -J x Tcste)] 
1+(Tcnocr—Tanocr)(Gs/Gs.nocr)U x Imp ste/T&) 


cS 


(18) 

3.3.1.2.2. PV output power. In this study, we will use one of the 

simplest models to predict the PV array output power, which 
assumes that: 


G 
Ppy = Ypvf py =) [1 +ap(Tc— Tcstc) | (19) 
T.stc 


where Ypy, fpy, Grstc, @p, Tc and Tcstc are respectively the rated 
capacity of the PV array, meaning its power output under standard 
test conditions, the PV derating factor (it is used to account for 
such factors as shading, snow cover, aging, and so forth), the 
incident radiation at standard test conditions, the temperature 
coefficient of power (which indicates how strongly the PV array 
power output depends on the cell temperature), the PV cell 
temperature in the current time step and finally the PV cell 
temperature under standard test conditions. 


3.3.2. Wind model 

3.3.2.1. Weibull distribution. The Weibull distribution is often used 
to characterize wind regimes because it has been found to provide 
a good fit with measured wind data [57-61]. The probability 
density function is given as follows [62-69]: 


K/V! K 
ep [—(V/c)"] 
fV)= cle) e (20) 


where V, K and C are respectively the wind speed, the shape factor 
(unitless) and the scale parameter (m/s). 

There are several methods which have been proposed to 
estimate Weibull parameters [70-72]; graphic method, maximum 
likelihood method and moment method are commonly used to 
estimate Weibull parameters. Also, a recent method proposed by 
Akdag and Dinler [73] in 2009 that is referred to “the power 
density method” will be exhibited in the following section. 

Akdag and Dinler [73] have suggested a new method, namely 
power density method, which is very useful to estimate scale and 
shape parameters thanks to its easy completion, simple formulation 
and also needs less computation. Indeed, it does not require binning 
and solving linear least square problem or iterative procedure. 

According to the PD method, the shape parameter, K, can be 
calculated using the following expression: 


(21) 


where Ep is the energy pattern which according to the literature is 
between 1.45 and 4.4 for overall wind distribution in the world [73]. 
The energy pattern can be expressed also as follows: 


ve 

v 
where (V)? and V? are respectively the average cube of wind speed 
and the average of wind speed cubes. 


After computing the shape parameter using (Eq. (21)), the scale 
parameter can be also calculated easily using the following expression: 


Ep = (22) 


V 


=PO) 23) 


where V, I” are respectively the mean wind speed and the gamma 
function. 


3.3.2.2. Energy and power in the wind. The kinetic energy, Ep, of a 
mass m of air that moves through a cross section A perpendicular 
to the wind speed V may be expressed as follows [74]: 


Pes jmv? (24) 


As the mass of air is the product of the air density p and the 
volume of the air that passes through the area A during the period t. 
The energy Ek may, therefore, be expressed as follows: 


Bye MpAVV? = = pAV? t (25) 
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It is well known that the power of the wind that flows at speed 
V through a blade swept area A increases as the cubic of its velocity 
and is given by [75-77]: 


P(V) = Av? (26) 


The wind power density, P4(p, V), defined as the power per unit of 
area perpendicular to the direction from which the wind is 
blowing. Pg(p,V) depends on the air density and wind speed. 
It is given as follows [78]: 


1 
Pa(p.V)=5P Vv? (27) 


In the above expression, if p is expressed in kg/m? and V is 
expressed in m/s Pg(p,V) is obtained in W/m?. P4(p, V)q), is the 
basic unit for measuring the power contained in the wind [79]. 

Wind power density of a site based on a Weibull probability 
density function can be expressed as follows [80]: 


Poe -i pikt 
a= [ PUY) AV = Spc re) (28) 


Once wind power density of a site is given, the wind energy 
density for a desired duration (a month or a year) can be expressed as: 
E 1. 3,/k+3 
ie alae ea )r (29) 


where T is the time period (or duration). For example, T is 720 h for 
monthly duration. 


3.3.2.3. Wind output energy. Basing on Weibull representative 
data, the output energy, Ewrr , is given by [81]: 


N 
Ewrr= © Pouti x f(Vi) x N x At; (30) 
iZi 


where Pouti, f(Vi), N and At; are respectively the fitted power 
output of a given wind turbine corresponding to the V; wind 
speed, the Weibull frequency distribution of V; wind speed, the 
total number of measurement for a particular period of the year 
and the hourly time interval. 


3.3.3. Diesel generator model 

The use of diesel generator is common in many hybrid energy systems 
to guarantee supply continuity and battery charge for energy storage. 
Experimental studies [21] have been found that for a diesel generator, 
a linear function fits for light load working conditions, while in 
proximity of rated power, it asks for quadratic expression. For a given 
time interval, the rate of fuel F, consumed by the diesel generator 
delivering the power P, is expressed in general as: 


F=aP*+bP+c (31) 


Table 1 
Design variables used for the studied a hybrid energy system. 


Variable Value 

Annual interest rate 6% 

Life span of the system 20 years 

PV panel cost (Kyocera KD 185 GX-LP) €437 

Solar panel installation fee 50% of the price 
Wind turbine cost (Whisper 200) €4379 

Wind turbine installation fee 25% of the price 
Unit cost of battery (Enersol 50) €142 

Backup generator cost (BC Electronic) € 900 

Usage% of battery's rated capacity 80% 

Battery’s rated capacity 624 W/h 
Battery’s life span 1500 cycles 


Maintenance cost for PV array 
Maintenance cost for wind turbine 


€0.005 per kWh 
€0.02 per kWh 


Table 2 


Hourly electrical consumption distribution for the different household appliances: case A. 
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where a, b and c are the coefficients of diesel generator as obtained 
from the manufacturers’ data. 


3.3.4. Battery model 

In hybrid energy systems, battery bank is required for storage 
of intermittent wind and solar energy. The battery state charge can 
be computed as follows: 

Battery charging, 


P,(t) = P,(t—1) x 1-0) 4 (20-20) xMp (32) 


Battery discharging, 


P,(t) = Pp(t—1) x (1-9) ee PO) (33) 
I 

where P;(t—1), P,(t), are respectively the battery energy at the 
beginning and the end of interval t, P;(t) and P;,(t) are respectively the 
load demand and the hybrid energy generated by PV array and wind 
generators at the time t, ø, n; and 77, are respectively the self-discharge 
factor, the inverter efficiency and the battery charging which are 
specified by the manufacturers’ data. 

The number of batteries Ngar can be computed using the 
following expression [82]: 


Nepat(Npv, Nwr) = Roundup | = | (34) 

Spat 
where Roundup(.), Sreq and Spar are respectively the function 
which returns a number rounded up to an integer number, the 
required storage capacity; the usage percent of rated capacity 
which guarantees battery’s life span and the rated capacity of each 
battery. We note that the required storage capacity Speq iS a 
function of number of solar panels Npy and number of wind 


3500 


turbines Nwr as follows: 


Max t f t Min t t t t 
Srec(Nev.Nwr)= È (Phy +Pwr—PpomdD At- È (Ppy +Pwr—Ppma)At 
t=1 t=1 
(35) 


where Maxt is the time when cumulative energy (kW h) is high- 
est; Mint is the time when cumulative energy (kW h) is lowest; At 
is the unit time step; Phy is the output power (kW) generated by 
the PV panels at time t; Pi; is the power (kW) generated by wind 
turbines at time t; and P5,,4 is the power (kW) demanded at time t. 
Besides, the PV power generation can be calculated as follows: 


Poy = Nev x P py Each (36) 


where P py gacn is the power (kW) generated by each PV panel at 
time t. P py gach Can be obtained as it was mentioned in Section 
3.3.1.2.1. Also, the power generated by wind turbines can be 
calculated as follows: 


Pwr =Nwr x P wr Each (37) 


where P wr _gach is the power (kW) generated by each wind turbine 
at time t. P wr cach Can be calculated using Section 3.3.2.3. 


4. Optimization formulation 


The objective function of the hybrid energy system design is 
the minimization of total design cost Cr which consists of total 
capital cost Cept and total maintenance cost Cynr as follows: 


Minimize Cr = Cepe+Cunr (38) 


We note that while the capital cost occurs in the beginning of a 
project, the maintenance cost occurs along the project life. 
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Fig. 5. Hourly profile of the household electricity load (case A). 
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Table 3 


Hourly electrical consumption distribution for the different household appliances: case B. 
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Consequently, costs at different times cannot be directly com- 
pared, but should first be made equivalent through the use of 
discount factors that convert a monetary value at one time to an 
equivalent value at another time [83,84]. For this reason, the initial 
capital cost 
P is converted into annual capital cost A using the following 
capital-recovery factor: 
A id+i" 

P (14i"+1 

where i and n are respectively the annual interest rate and the life 
span of the system (in years). Therefore, total annual capital cost 
Ccpr can be given as follows: 


(39) 


A 
Ccpt = P [NpevCpy +NwrCwr + NpgarCear + Cpackup] (40) 


where Npy is number of PV panels, Cpy is unit cost of PV panel, 
Nwr is number of wind turbines, Cwr is unit cost of wind turbine, 
Ngat is number of batteries, Cgar is unit cost of battery; and 
Cgackup İS cost of backup generator for the use when the electricity 
demand is greater that the PV, wind and stored energies. The unit 
cost of solar panel Cpy consists of panel price and installation fee; 
and the unit cost of wind turbine Cwr consists of turbine price and 
installation fee. 

The total annual maintenance cost Cynr can be expressed as 
the following: 


_ [ppv 24 ot wT 4 pot 
Cunt = |Cunr x È (Ppy x At)+Cynr x È (Pwr x At)| x 365 (41) 
t=1 t=1 


where Cr and Cyiny are respectively the maintenance cost per 
kW h for PV array and wind turbine. Here, we have to note that 
due to the vulnerability of battery in the renewable power 
generation system, the replacement cost of the battery may be 
also included in the objective function of the total design cost. 

In the optimization problem, several constraints may be 
imposed. First, the total energy amount generated by the PV 
panels and wind turbines should be greater than or equal to total 
energy amount required by users as follows: 

24 24 24 
È (Pov x At)+ X (Pwr x At)= X (Poma x At) (42) 
t=1 t=1 t=1 

Second, the numbers of PV panels and wind turbines should be 
non-negative integer variables as the following: 


Npy = Integer, Npy >O (43) 


Nwr=Integer, Nrw>0 (44) 


Third, the problem can have optional constraints of resource 
limitation as follows: 


Npy < Npa“ (45) 
Nwr < Wr (46) 
Neat < Near (47) 


where Npy*, Nw and Ngay are respectively the maximum avail- 
able number of PV panel; the maximum available number of wind 
turbines and the maximum available number of storage batteries. 

Table 1 provides the values of the used design variables. It is 
noticeable that because the life span of each battery is 4 years, five 
times of Ngar is required to satisfy the life span of the PV-wind system. 

In this paper, a new proposed method was used to solve the 
optimization problem. This method has a very smart procedure 
exploration of the parameter space. In fact, this method is 
relatively simple when comparing with dynamic simulation. 
It explores all the variable space and records the best element. 
The objective function Cr is evaluated at a point under the 


S. El Alimi et al. / Renewable and Sustainable Energy Reviews 37 (2014) 408-423 417 


considered constraints. The new value of Cr is compared to the 
previous one. If it is better than the previous one, this value is 
stored as well as the corresponding solution, then the process 
continues. Otherwise, the process retains the previous point and 
the procedure continues until stopping conditions are reached. 
The proposed method draws the region where the solution is 
acceptable in the parameter space. The flowchart of ascertaining 
the optimal system configuration in our case study is as the 


following: 


Get the input load data, solar radiation, 
wind velocity and also the techno 


economic data of HSE components 


Ascertain hourly power generated by single units of 
renewable generators using mentioned models 


Npy = Npy +1 


If No 


——— |" yes 


Calculate Sreq, Near, Cr 


Cr ont = CTNpvopt = Noy 
NwTopt = Nwr 
— 
If No 
If yes 


Flowchart of ascertaining the optimal system configuration. 


where Ap=(Ppy+Pwr)—Ppma is the difference between the 
system output power and the power load consumption for each 
time step; 


© Crop: optimal actualized total cost of configuration system; 
© Npvopt: optimal number of photovoltaic modules; 
© Nwropt: optimal number of wind turbines. 


The above mentioned optimization technique use mathemati- 
cal methods to model the PV, wind and battery sources. At the 
same time, this model depends on availability of weather informa- 
tion and details of manufacturer’s data. Indeed, the optimization 
model requires the following inputs: 


— Resources availability (solar radiation data, clearness index and 
wind velocity) of the site in question for each computational 
time step; 

— Electrical load distribution; 

— Technological and economic inputs: all characteristics and 
prices of components. 

— PV module: orientation and tilt angles, optimal output power, 
thermal power coefficient, temperature of the cells functioning, 
aging coefficient, module unit cost, installation fee, mainte- 
nance cost per kW h and life span. 

— Wind turbines: output power curve, hub height, wind turbine unit 
cost, installation fee, maintenance cost per kW h and life span. 

— Battery: nominal capacity, round trip efficiency, float life, unit cost. 

— Backup generator: empirical consumption rate coefficients, unit 
cost, life time. 


The proposed optimization algorithm simplifies the task of 
evaluating system design for variety of combinations. This model 
simulates the operating of system by making energy balance 
computations in each time step. In fact, for each time step, this 
model compares the electrical demand and the output energy 
that system can supply in that time step, and calculate the flows 
of energy to and from each component of system under the 
considered constraints. It then runs for each possible combina- 
tion of the system (configuration). For configurations that include 
batteries and backup generators, we also evaluate in each time 
step how to operate generators and calculate the required storage 
capacity. The energy calculations for each time system config- 
uration determines whether it is a feasible one; means that 
whether it can meet the electric demand under the conditions 
and the constraints that we specify, and it estimates the actua- 
lized annual cost system over its lifetime. This system cost 
calculations account for costs such as capital, replacement, 
maintenance, installation and interest. After simulating all fea- 
sible configurations, only the best annualized cost-effectiveness, 
which minimizes the objective function C;, will be stored. We 
note that we display all optimal configurations of each possible 
hybrid combination. 


5. Outcome results 


This section covers the three load consumption distributions, 
the output power computations of the available resources and the 
techno-economic deigns of the optimal HES configurations. 


5.1. Load consumption 


In this study, three cases studies will be considered to define 
the power consumption of an urban residential building located in 
the region of Monastir, Tunisia. The latitude and longitude of this 
town are respectively 36°N and 12°E and an altitude of 10 m above 
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sea level. We note that this city has a good wind resource and solar 
deposit that promote the use of renewable energy resources for 
electricity generation. 

For the load consumption, a total of 22,054 W/h is considered 
as the average daily energy demand for a typical house in the city 
of Monastir. This average value of a daily energy load was 
confirmed by a statistical survey which has been established by 
the Electricity and Gas Tunisian Company (STEG). 

Three different distributions are defined, the first one (case A) 
corresponds to the families who leave their habitats from 8 am, a 
half-round home at mid-day, another out of the house at 2 pm and 
the whole family back home at 6 pm. The peak consumption has 
the values of 416.1, 3019.1 3261.1 W which occurred at 8 pm, noon 
and 8 am, respectively. Table 2 and Fig. 5 show, respectively, the 
distribution of the household electrical charge for the different 
energies equipments and its corresponding hourly profile. 

The second profile distribution of the house electricity demand 
(case B) corresponds to the families who leave their habitats from 
8 am and return back at 2 pm. This path can explain the first peak 
consumption at 6 am with the value of 619 W and during the 
second half of the day and exactly at 3 pm and 4 pm it reaches a 
consumption of 2916 W. The last peak of the electricity demand 
take place during the hours between 8 pm and 10 pm with a 
maximum of 3261 W. The energy consumption often occurs in the 
afternoon especially for the needs of air conditioning and in the 
night for almost imperative lighting and operation of different 
appliances. Table 3 and Fig. 6 show respectively the distribution of 
the household electrical charge of the different energy equipments 
and its corresponding hourly profile. 

The third distribution of the power consumption (case C) 
corresponds to the families which are using for almost of the time 
their domestic electrical appliances at the night. For this reason, 
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after slight morning peak consumption reached at 6 am having the 
value of 619 W, the electric charge increases from 6 pm achieving 
another peak load corresponding to 4261 W at 8 pm. Table 4 and 
Fig. 7 show respectively the distribution of the household elec- 
trical charge of the different energy equipments and its corre- 
sponding hourly profile (Table 5). 


5.2. Output power resources 


Information about the average hourly annual profile of the 
wind speed and data interpolation of the wind power curve 
(Whisper 200 W) are very necessary for a confident estimation 
of the wind turbines efficiency (Fig. 8). 

Using the average hourly annual wind speed profile corre- 
sponding to the studied region, the interpolated wind power curve 
of Whisper 200 W and the Weibull model for estimating the wind 
power (Eq. (15)), we obtain the average hourly annual wind 
production distribution which is shown in Fig. 9. 

Using the HDKR model for estimating global solar radiation and 
the PV conversion model (Eq. (7)), the average hourly annual PV 
output power can be shown in Fig. 10. 

We note that this study aims mainly to investigate the load 
profile distribution impact on the optimal hybrid energy config- 
uration. The use of an annual average of solar radiation and wind 
speed is valid in the North of Africa and especially in Tunisia 
because this it has quite stable climatic conditions. However, for 
other locations, especially for higher latitudes where there are 
extremely high seasonal variation, especially for higher latitudes, 
it would be interesting to apply the proposed method using an 
average hourly meteorological data. 
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Fig. 6. Hourly profile of the household electricity load (case B). 


Table 4 


Hourly electrical consumption distribution for the different household appliances: case C. 
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5.3. Optimal HES configurations 


In this simulation, for each case of residential electricity 
demand set (A, B or C) seven configurations are possible (2"~! 
where n is the number of available sources in the SEH (PV/ 
Wind/Diesel). For each hour, the developed code simulates all 
possible combinations. Hour after hour, the load consumption is 
compared to the available potential energy and the energy flow 
is adjusted for each system component. Besides, for each hour, 
we can monitor the functioning of the diesel generator, and the 
charging and discharging of the batteries bank. As the hybrid 
system model is introduced with the technological options 
availability, component costs and resources availability, these 
inputs are simulated for different combinations of possible 
components. Hence, Energy balance calculations are performed 
for each configuration of hybrid system. Profitability operation 
is assessed on the life of the project which is estimated to 20 
years. The used optimization technique is very interesting for 
combinatorial systems’ optimization. Indeed, the technique uses 
an enumeration of all candidate solutions process ensuring the 
provision of a global optimum. This method takes the following 
objective function: 


ERS A 
Minimize Cr = P [Npv Cpy +NwrCwr+NepatCpat + Cpackup] 


+ [Chin x 24 (Poy x Ab+C\nt x 524 (Pwr x An] x 365 


The blown table shows the optimal hybrid configuration for 
each case of distribution for the residential electricity consump- 
tion in the region of Monastir. The total annual cost and the 
levelized cost of energy of each configuration of the hybrid 
power system are recapitulated in this table. Here, we note 
that the levelized cost of energy is the constant unit cost 
(per kW h) of the capital annualized cost of the system over 
its life. It is clear that for the same average daily energy 
requirement load, the optimal hybrid energy system is not the 
same because of the difference in the distribution of peak 
consumption during the simulation of the three cases of resi- 
dential load profile. 

We note that only the first case (Case A) PV power generation 
can contribute in the global energy production due to the two 
electrical peak demands which occur at 6 am and 2 pm. In this 
case, the presence of 7 PV modules is very effective for both 
their maximum production coincides with the highest electrical 
demand. We note also that this configuration is the most cost- 
effective one with a levelized cost of energy equal to 0.3082 € 
per kW h over 20 years. Then, we can deduce that the combina- 
tion (PV/wind) is the best energizing solution because it satis- 
fies, with the minimum annual cost (2481.1023 €), the worst 
case power distribution in terms of stability. In addition, the 
optimal configuration of the hybrid energy system in the second 
case (Case B) does not longer contain the PV solution. This is 
explained by the fact that periods of high electricity consump- 
tion (3-4 pm and 8-10 pm) does not correspond to the 
maximum PV output power (see Fig. 10), which can lead to 
quite high operating generators cycles. Hence, the optimal 
configuration is composed by three wind turbines which their 
maximum output power fulfill at 3 pm and 4pm (see Fig. 9) and 
15 storage batteries. This second optimal configuration gener- 
ates a cash surplus because the levelized cost of energy is now 
evaluated to 0.3101€ per kWh. Finally, the third optimal 
configuration of HES for the third electric charge distribution 
(case C) retains the same optimal wind design but the number of 
batteries storage has been increased by 8 ones in comparison 
relatively with the second case (Case B). This increase of storage 
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Fig. 7. Hourly profile of the household electricity load (case C). 


Table 5 
Techno-economic design of each optimal HES configuration. 


Load case Optimal Annualized total Levelized cost 
configuration cost Cr (€) of energy (€/kW h) 

Case A 2 WT-7 PVM-15 BS 2481.1023 0.3082 

Case B 3 WT-0 PVM-15 BS 2496.3814 0.3101 

Case C 3 WT-0 PVM-23 BS 2991.5896 0.3716 


E.I.: WT: Wind Turbine, PVM: Photovoltaic Module and BS: Battery Storage. 


capacity is explained by the dominant concentration of electri- 
city demand during the night hours that's why the resulted total 
annual cost and the levelized cost of energy become higher than 
the last two cases (which are respectively 0.3716 € per kW h and 
2991.5896 €). 

The main conclusions drawn from this investigation can be the 
following: 


© The electricity demand distribution is an important factor in 
the hybrid renewable energy system’s optimization and not 
only the resources availability and the quality of the local 
renewable deposits. 

e@ The developed optimization model has proved its efficiency in 
terms of computation time and reliability to achieve the overall 
optimal configuration without the risk of falling on potential 
local solutions. 

© The HES (PV/wind/diesel) combined with electrochemical sto- 
rage banks can be a cost-effective energy system and econom- 
ical solution in medium and long term. The following table 
recapitulates the optimal hybrid configuration and its total 


annual cost and the levelized cost of energy of each load 
demand category. 


6. Conclusion 


In this study, the sizing of a (PV/Wind/Diesel) hybrid power 
system with electrochemical storage banks has been optimized 
using a proposed generic technique based on the principle of the 
determinist methods. The main objective is to minimize the annual 
operating cost and the produced Kilowatt of the system over 20 
years. First, a complete HES modeling has been developed. Second, 
the modeling of the supplied electric auxiliary power (diesel gen- 
erator) and its fuel consumption with the instantaneous state of 
charge and discharge battery banks has been implemented. 

The effectiveness analysis of the SEH was performed in order to 
decide which investments will be beneficial to cover the energy 
consumption of standard residential buildings taking into account 
the non-uniform distribution of the electric charge in the city of 
Monastir, Tunisia. 

An optimization model of the actualized annual cost over the 
HES life span was made. This optimization model is very useful for 
the simulation of the medium and the large scale of installed 
capacity for autonomous renewable systems for. Indeed, this model 
starts with the optimization design from the optimal energy 
resources evaluation (estimated solar flux correlations, optimal tilt 
and azimuth panels, optimal wind speed distribution, optimal 
interpolated power curves...and so forth) until the development of 
real relationships between several input variables of the system 
(number of PV modules, number of wind turbines and storage 
batteries) and output variables. Then, the results showed that the 
optimum hybrid system for the same electricity demand depends on 
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Fig. 10. Hourly output power produced per unit of PV module (Kyocera 185 GK-LP). 


the probability distribution of the daily load consumption. The 
combination (PV/Wind) is the best way to satisfy the most instable 
electrical demand distribution and in this case the levelized cost of 
energy does not exceed the value of 0.3082 € per kW h. For electrical 
applications for which the peak loads occurs after noon or night, 
wind energy is the promising renewable resource by excellence 
under the local climatic conditions of the studied region. In fact, the 
optimal levelized cost of energy is ranging from 0.3101 and 0.3716 € 
per kWh. For these cases, the load distribution profile and the 
storage capacity will make the difference. 

Finally, the optimization methodology described in this study 
provides a very important systematic approach to the design, 
analysis and optimization of hybrid energy systems. 
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